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Abstract 

Static magnetic quadrupole moment of a nucleus, induced by T- and P-odd nucleon- 
nucleon interaction, is investigated in the single-particle approximation. Models are 
considered allowing for analytical solution. The problem is also treated numerically 
in a Woods-Saxon potential with spin-orbit interaction. The stability of results is 
discussed. 
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1. Magnetic quadrupole moment is a static characteristic of a quantum system which is 
forbidden by P- and T-invariance. Nuclear magnetic quadrupole moment (NMQM) can be 
induced both by the nucleon electric dipole moment]]]] and by P- and T-odd nuclear forces 0. 
The interest to NMQM is due to the experimental searches for P- and T-odd effects in atoms 
and molecules (see, e.g., book||). 

The manifestation of electric dipole moment (EDM), which also violates P- and T- 
invariance, in atomic and molecular phenomena is strongly hampered by the electrostatic 
screening. In a stationary system of nonrelativistic pointlike particles interacting via Coulomb 
forces such a screening of average electric field acting on any paricle is complete. Therefore, 
in such a system for a particle EDM there is nothing to interact with, which means that this 
EDM just cannot be observed @, §]. 

The nuclear dipole moment becomes observable however due to the finite size of a nucleus, 
more exactly due to different distribution of its charge and EDM||. One more way to 
transfer nuclear P- and T- violation to an atom or molecule is via NMQM0. It demands of 
course nuclei with spin I > 1/2. Besides, the NMQM induces P- and T-odd effects only in 
atoms (molecules) with unpaired electron angular momenta since it interacts directly with 
magnetic field of the electrons. However, when operative, the NMQM is much more effective 
for circumventing the electrostatic shielding in atoms and molecules |], [2|, |3|]. 

It has been shown in Ref.0 that the NMQM induced by the P- and T-odd internu- 
cleon interaction can be much larger than that due to the nucleon EDM. In that paper the 
quadrupole moments, induced by that interaction, were evaluated in a simple model where 
the profile of nuclear density was assumed to coincide with that of nuclear potential. In the 
present article we calculate NMQM within a more accurate approach. Namely, we use a 
realistic description of the nuclear density; the nucleon wave functions and Green's functions 
are obtained with Woods-Saxon potential which includes the spin-orbit interaction. We in- 
clude also the contribution of the current generated by the spin-orbit interaction; contrary 
to naive expectations, this contribution does exist for an outer neutron, but does not for 
an outer proton. One more model admitting a closed analytical solution is considered, that 
of the oscillator potential. We restrict throughout the present paper to the single-particle 
approximation, that of a valence nucleon above a spherically-symmetric core. 

This approach was recently used by us[[7| for treatment of nuclear anapole moments, 
P-odd, but T-even characteristic. 

2. Let us begin with discussing the T-and P-odd nucleon-nucleon potential. In the local 
limit and to first order in the nucleon velocities p/m it can be written as follows (see, e.g., 
book§) 

G 1 -> 

W 2 = ^—Y. ((Ufa ~ ibaVb) ■ V5(f a - r b ) 

H'ab [°a X &b] ■ {(Pa~ Pb),S(f a ~ f b )}) (1) 

where the notation { , } means anticommutator. The dimensionless constants £ char- 
acterize the magnitude of the interaction in units of the Fermi weak interaction constant 
G = lCT 5 /m 2 and are supplied with subscripts in order to distinguish between protons and 
neutrons. 
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After averaging this expression over the core nucleons we obtain the P- and T-odd mean 
field potential for an outer nucleon 

Here p(r) is the density of the core nucleons normalized by the condition / dfp(r) = A 

04 »i); 

Z N 

Sa S ap + San i 

the subscript a takes the values p and n for an outer proton and neutron, respectively. 

Let us note that, as distinct from the case of the P-odd, T-even interaction, no contact 
current is generated here in the single-particle approximation, even if one starts from the 
two-body interaction ([!]). Indeed, it is only the last term in ([]]), dependent on £' ab , which 
contributes to the contact current operator 

J C = ^E{[^,e a f a ],5(f-f a )}. (3) 

However, even this contribution vanishes obviously after averaging over the core nucleons. 

Now, the correction 5^/ to the valence nucleon wave function generated by the interaction 
(|) is a solution of the equation 

(H - E)8V(r) = -W(r)V(r), (4) 

where H an d \P(r) are the unperturbed mean field Hamiltonian and the unperturbed nucleon 
wave function. To begin with, let us discuss a simple model where the profiles of the nuclear 
density and the central mean field potential coincide, and the spin-orbit potential is absent 

p(r) = -£-U(r). 

Eq.(f|) transforms here as follows: 

(H - E)5V(r) = ^G=±.£±[H , a ■ ft* ft, (5) 

which gives 

Here Qu m is a spherical spinor, R(r) is the unperturbed radial wave function of a nucleon, 
and K= (l-I)(2I+l). 

Even in a more general case, beyond this model, it is convenient to define the correction 
5R(r) to radial wave function by the following relation: 

**(r) = -^po(B-m Ilm {n)8R{r). (7) 
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The correction 5R(r) can be calculated using two independent solutions of the radial Schrodinger 
equation U\{r) and u 2 (r), regular at the origin and at the infinity, respectively. These solu- 
tions are normalized to the unit Wronskian: 

du2 du\ 

Mi — — u 2 = 1. 

dr dr 

This correction is 

6R ( r ) = — I dr u W'{r )~^r u o(r ) - — J o dr u 1>IV {r )-^rM r h ( 8 ) 

where V = 21 — I, u (r) = rR(r) and /(r) is a density profile /(r) = p(r)/p . 

The magnetic quadrupole moment operator My is defined by analogy with the electric 
quadrupole one Qij, via the interaction with the corresponding field gradient: 

Hq = —QijViE,, (9) 
H M = —-MijViBj. 

The symmetric tensor M^- is related in the following way to the current density J n : 

M^ij J dr(ri6j mn ~\~ r j€i mn )r m J n . (10) 

For a valence nucleon this operator can be presented as|2|. [3|] 

M i:j = (zpiriVj + TjOi - -Sij(a ■ f)) + 2q(r i l j + rjl^ , (11) 

where p is the nucleon magnetic moment, and q is equal 1 for a proton and for neutron. 
With the usual definition 

M = (Im = I\M zz \Im = I) 
one obtains after taking expectation value over angular variables 

M=^ Po ^-(p- q )*L^l(5R\r\R). (12) 

The radial matrix element here is 

POO 

(5R\r\R)= r 2 dr8R{r)rR(r). (13) 
Jo 

For the simple model described above (§) the matrix element can be calculated analytically 
with the following result 

< M ' r '«) = (14) 
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One more model allowing for an exact analytical result for MQM is that of the oscilla- 
tor potential. Here it is convenient to start from expression ([13]) for the matrix element. 
Separating the tensor structure in (reftvpot), QTTT) one obtains 

{5Rm = _ J_ V (Olf • n)\n)(n\a ■ f|0) + (Ojg ■ gjf^HQg ; n)|0) ^ 
4m „ E — E n 

For a harmonic oscillator 

(^■rO = - L a[(^-^),fll- 

mur 

Substituting this identity into (|I5D and using the completeness relation, we find 

(5R\r\R) = -^- 2 (0\[f(a ■ ft), a ■ p]|0). 
Taking again expectation value over angular variables, we obtain 

3. Expression ( fTT| ) for the MQM corresponds to the contribution of the convection and 
the spin electromagnetic current densities. In this section we discuss one more contribution 
to MQM, that originating from the momentum dependence of the spin-orbit two-nucleon 
forces. In the single-particle approximation this current density is0 

j T pn Ndf(r) _ 
3fs^eU^p —-^-axn (17) 



for a valence proton and 



It = -eU? s n j Po f(r) V x (^(^(f)) 



for a valence neutron. Here UfJ 1 is the constant entering two-body proton-neutron spin-orbit 
treated in the contact limit: 

Uu = \ E U l n (Pp - Pn) -(a P + a n )x V5{f p - f n ). (19) 

^ pn 

The proton-proton spin-orbit interaction does not contribute to the current density in the 
zero-range limit we use. 

The direct calculation shows that, contrary to naive expectations, the proton spin-orbit 
current fllTD does not contribute to the static NMQM. However, for an outer neutron the 
corresponding correction does not vanish. It equals 

M£ = -2^ZdeU?jj^±(8R\rf{r)\R). (20) 

4. We are ready now for a more realistic single-particle calculation. This numerical 
treatment is based on the Woods-Saxon potential including spin-orbit interaction and on a 
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realistic description of nuclear density. The profiles of both density and the central part of 
nuclear potential are described by a Fermi-type function 

/M = i + eM^y < 21 > 

The total single-particle potential U (r) is chosen in a standard Woods-Saxon form 

U{r) = U f(r) + U ls -^P-(fa) + U a (r), (22) 
r ar 

where Uc{r) is the Coulomb potential of a uniformly charged sphere. 
We use the values of the density parameters from book||: 

R = 1.11 A 1/3 fm; a = 0.54 fm p = 0.17 fm~ 3 . (23) 

The Woods-Saxon potential is parametrized as in@: 

R = R ls = 1.24A 1/3 fm, a = a ls = 0.63 fm, 

N - Z N - Z 

U = (-53.3 ± 33.6 — - — ) MeV, U u = -0.263 (1 + 2 — — ) U (24) 



The spin-orbit interaction constant as fitted in Ref.[lC 



is: 

U[ s n = 134.3 MeV ■ fm 5 . (25) 

The correction 5R(r) calculated in this way is plotted in Fig.l together with the model 
function @. Obviously, the latter is a reasonably good approximation to the correction 5R, 
as calculated numerically in the more realistic approach. 

Our results are conveniently presented in terms of the dimensionless constant r related 
to the NMQM as follows 

This constant itself consists in general of two contributions: 

r = r + T ls , 

r = ^-(fM-q)2(5R\r\R), (27) 

m p 

rll = -4U^pl(5R\rf(r)\R), (28) 



A' 

0. 
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The values of r calculated for two neighbouring nuclei, with odd Z and odd N respectively, 
are presented in Table 1. 

5. In conclusion let us compare the results obtained for the NMQM, generated by P- and 
T-odd potential, with those for the nuclear anapole moment (AM), generated by P- violating, 
but T-even potential. In particular, we wish to compare the stability of nuclear single-particle 
calculations for those two moments, T-odd and T-even. 

For the sake of comparison with the constant r calculated here, it is convenient to delete 
from the dimensionless AM characteristic k the fine structure constant a (related to the 
electromagnetic interaction of an atomic electron with nuclear AM) and the P-odd nucleon- 
nucleon constant g (the T-even analogue of the constant £ used here). The typical value of 
this AM characteristic is JTT], |3], |7j 

J^~J^A 2 '\ (29) 
ag mr 

where r = 1.2 fm. 

As to r, its typical value is 

r ~ fiK4m 2 U rl (30) 
The ratio of those two factors, i.e., of the T-odd effect to T-even one, 

K 9 A~ 2/3 ~ 0.15 K A' 2 '* (31) 
4m£7o?o 

is very small. The origin of the AM enhancement ~ A 2 ^ can be traced back to the AM 
dependence on the geometrical cross-section of nucleus, this is a bulk nucleus effect [[□], ||, [7J. 
As to the NMQM, its magnitude depends completely on the nuclear boundary (see eqs.(H), 
©)■ 

It results not only in the relative suppression of the T-odd effect. The value of NMQM is 
more sensitive to the details of the nuclear model than that of AM, it is less stable. 

However reliable theoretical predictions both for AM and NMQM can be obtained only 
when the single-particle calculations will be supplemented by a serious treatment of many- 
body effects. 

This investigation was financially supported by the Russian Foundation for Fundamental 
Research. 
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Nucleus 



—Upo/Uo Harmonic oscillator Woods-Saxon 



133 Cs (lg p 7/2 ) 


T = T 


0.16 


0.26 


0.18 


137 Ba (2d- /2 ) 


To 


-0.09 


-0.17 


-0.12 




Tls 


-0.02 




-0.02 




T 


-0.11 




-0.14 



Table 1: The dimensionless MQM as calculated in different approaches 

Figure Caption 

Figure 1. The correction 5R(r) for 137 Cs. Dashed line is the model function (|6]). Full line 
is the SR in Woods-Saxon potential 
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